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ABSTRACT 

We present high angular and spectral resolution HI 21 cm line observations toward the cometary- 
shaped compact HII region G213. 880-11. 837 in the GGD 14 complex. The kinematics and morphology 
of the photodissociated region, traced by the HI line emission, reveal that the neutral gas is part of an 
expanding flow. The kinematics of the HI gas along the major axis of G213. 880-11. 837 shows that the 
emission is very extended toward the SE direction, reaching LSR radial velocities in the tail of about 
14 km s~^. The ambient LSR radial velocity of the molecular gas is 11.5 km s~^, which suggests 
a champagne flow of the HI gas. This is the second (after Gill. 61+0. 37) cometary HII/HI region 
known. 

Subject headings: HII regions, ISM: individual objects (GGD 12-15, GGD 14, G213. 880-11.837), radio 
lines: ISM 



1. INTRODUCTION 

The spatial morphology, abundance, and kinemat- 
ics of the neutral gas around compact HII (CHII) 
regions and its relation with other physical pa- 
rameters remains poorly studied for most HII re- 
gions. To date, only a few photodissociated re- 
gions (PDR) have been observed at radio w avelengths 
with high angu lar resolutio n (e.g. [l/ebron fc Rod riguezj 
1997tlG6me"zet al. 1998; Ga rav et al.lll998 : Broga n et al.l 
1999t iLebronet al.. i2001i: iCappa et al.l 12008. ). GGD 
14 is an active star forming region, part of the 
group of red nebulous objects GGD 12/ 13/14/15 
(|Gvulbudaghian. Glushkov. fc DenisvukI I1978D embed- 
ded in th e Monoceros molecula r cloud at a distance of 
d^l kpc ([Rodriguez et al.|[l980[ ). There is an extended 
(^1.8 pc) bipolar outflow in the region, traced by CO 
emiss ion with a m ajor axi s along the NW - SE direc- 
tion jRodriguez et al. 1982:: iLittle et all 119901: iQin et al.l 
I2008D . This bipolar outflow seems to be excited by 
a small region of free-free emissio n (named VLA 7 by 
iGoniez. Rodriguez fc Garavl |2000[ ). possibly a thermal 
iet (jEisloffel et al.ll2000t) . The powering source of this 
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bipolar outflow is also associated with an H2O maser 
(iRodriguez et al.lll982D. a hvdrogen molecular emission 
object (JFang &: Yaol 1200^ and coincident in position 
with a class I/O source, recently detected in the mid- 
infrared (IRS 9Mc: Sato et al.: 2008 ). 

The GGD 14 star forming region contains a clus- 
ter of radio sources dominated by the compact (^0.01 
pc) cometary-shaped HII region (VLA 1; G213.880- 
11.837), which is excited by a BO. 5 ZAMS star 
(IRodriguez e^.l 119801: iKurtz et all 119941: iGomez et aLJ 
[l998t iGomer^ drigue z fc Garavl l2000. 2002) and is as- 
sociated with IRAS 06084-0611. A detailed kinematic 
study of t he gas within the c ometary HII region, pre- 
sented by 'Go mez et al.l (|1998D , shows that the ionized 
gas is undergoing a champagne flow with the head of 
the ionized champagne flow located to the NW and the 
extended emission toward the SE. There is a velocity gra- 
dient along the major axis of the ionized champagne flow 
with the LSR radial velocity increasing from ~11 km s~^ 
at the head of the flow to ~ 15 km s~^ at the tail. In 
a champagne flow the velocity at the head is expected 
to be very similar to the velocity of the ambient molec- 
ular gas, and this is the case for the champagne flow in 
GGD 1 4 where the ambient velocity is about 11 km s~^ 
(|Torrel les et al. 1989; Oin et al. 2008). 

fe4^1a «^i*yJW'5^fcnaa^x|^^3izaA9i^ ^O'^' 

G6mez et at.. (,19981 1 reported the existence of an 
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unresolved PDR surrounding the ionized gas. The 
distribution of the molecular gas in the GGD 14 region 
is complex with the presence of differ ent veloci ty compo- 
nents ( Rodriguez et al. 1980; Giiste n fc Marca ide 1986; 
Torrelles et all 119891 ILittle et all 119901 : lAng ada et al.l 
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19961). Recent ^300(2^1) observations by 
1^08) toward the GGD 14 region show that the 
molecular emission peaks at a velocity of ~11.5 km s~^, 
that could be considered as the velocity of the ambient 
molecular gas. 

In the present work we analyze new HI 21 cm line 
observations toward the cometary HII region G213.880- 
11.837 in GGD 14 with higher angular (15") and spec- 
tral (~1 km s ~^) re solution than the data presented by 
iGomez et al.l ()1998l) . These new observations resolve the 
morphology of the HI region that surrounds the ionized 
material. In §2 we present the observations, in §3 the 
results, in §4 an interpretation of the HI observations, 
and in §5 the conclusions. 

2. OBSERVATIONS 

The HI 21 cm {v = 1420.406 MHz) observations to- 
ward the star forming region GGD 14 were made using 
the Very Large Array (VLA) of the NRAO^ The obser- 
vations were taken on December 23, 1998 in the C con- 
figuration (under project AL459). The total bandwidth 
of 0.781 MHz was centered at a Vlsr of 12 km s^^, with 
256 channels 0.64 km s~^ wide each (3.1 kHz). Hanning 
smoothing was applied to the line resulting in a spectral 
resolution of ^1.2 km s~^. The flux density scale was 
determined from observations of the amplitude calibra- 
tor 0134-f329, for which a flux density of 15.87 Jy was 
adopted. The phase calibrator was 0605-085 for which 
a bootstrapped 1.4 GHz flux density of 2.11 ± 0.01 Jy 
was obtained. The data were edited and calibrated fol- 
lowing standard procedures and images were made using 
the NRAO software AIPS. Line images were made by 
subtracting the continuum (line-free) channels from the 
visibility data using the task UVLSF. The images were 
made with the task IMAGR, using the ROBUST parame- 
ter equal to zero (Briggs 1995i) , resulting in a synthesized 
beam of 17 x 14 , P.A.=— 26°. The rms noise level in 
a single spectral line channel, after Hanning smoothing, 
is ~2.0 mJy beam~^. 

3. RESULTS 

The detected radio continuum emission at 21 cm, 
made with the line-free channels, is marginally resolved 
(4"3x2!'9, P.A.=148°), with a flux density of 72±2 mJy 
and peak position at a(2000)= 06^ 10™ 50f607 ±05004, 
(5(2000)=-06° 11' 50"10 ±0"05. The HI 21 cm line, is 
detected both in absorption and emission in the velocity 
range from '-~^1 to 2 km s~^, in agreem ent with previous 
results reported bv lGomez et al.l (|1998D . 

Figure 1 shows the HI 21 cm spectrum integrated over 
the whole sou rce in a box wi th size of ~l'3xl'3. As was 
noted by Go mez et al.l ()1998[ ). the HI absorption compo- 
nent is considerably broad towards blueshifted velocities, 
showing three minima and extending from velocities of 



-1 km s ^ to 



-12 km s . In Table 1 we present 



the results of a Gaussian fit for the three HI absorption 
components plus a single Gaussian for the emission com- 
ponent (the fit is shown as a dashed-line in Figure 1). In 
this Table we give the peak fiux (Si), the center velocity 
(Vls_r) and the FWHM (full width at half maximum) 
of the line (AV) for each Gaussian component. The HI 
absorption components are centered at Vlsr oi ^ 1, 6 
and 11 km s~^ (see Table 1). On the other hand, the 
HI emission is present in a velocity range from ~ 13 to 
20 km s~^, with the maximum centered at VLSi?=15.5 
km s~^ (Figure 1 and Table 1). 

Figure 2 shows individual channel images of the HI 21- 
cm line toward GGD 14. The unresolved HI absorption 
appears coincident with th e 1.4 GHz continuum peak 
posit ion (named VLA 1 by iGomez. Rodriguez fc Garavl 
|2000( ). marked with a cross. The triangle, towards the 
NE of the CHII region, indicates the radio continuum 
peak position of the radio source VLA 7, which has 
been proposed as the exciting source of the b ipolar out- 
flow (jGomez. Rodriguez fc Garavll2000L l2002[ ) and is co- 
incident with the H2O maser emission in the region 
(jRodriguez et al.lll980[ ). As observed in the HI line in 
emission, the PDR in GGD 14 appears spatially resolved 
with a deconvolved angular size of 45"x40", P.A.=34° 
('--^0.2 pc) and with a peak intensity spatially shifted to- 
ward the SE side of the cometary HII region (the PDR 
peaks at a(2000)= 06^^ 10™ 51f3, <5(2000)=-06° 11' 52" 
± 2", about 0.1 pc in projected separation from the peak 
of the CHII region) . 

4. DISCUSSION 
4.1. Physical parameters of the HI line 

The individual HI 21 cm channels, toward GGD 14 
(Figure 2), show an extended emission component to- 
ward the SE and a strong line absorption coincident with 
the cometary HII region. The PDR, traced by the HI re- 
gion, is thus being produced mainly by the BO. 5 ZAMS 
star that also ionizes the cometary HII region (VLAl; 
G213.880-11.837). 

It is possible to estimate the HI optical depth toward 
the CHII region using the HI absorption profile. As- 
suming that the continuum source is located behind a 
homogeneous HI cloud, t he observed line b rightness tem- 
perature is given by (see lRohlfs fc Wilsoru.2003. ): 



TL{v) = {fHiT,^-fofcT,){l-e- 



-rL{v)^ 



^ The National Radio Astronomy Observatory is operated by 
Associated Universities, Inc., under cooperative agreement with 
the National Science Foundation. 



where T^x is the excitation temperature of the HI, Tc 
is the continuum brightness temperature of the CHII re- 
gion at the frequency of the HI line, ini and fc are the 
beam filling factors of the HI and the continuum, fo is the 
fraction of the continuum source covered by the HI cloud 
and tl {v) is the optical depth of the HI gas at a velocity 
V. In this case, we consider that the continuum source 
is completely covered by the HI cloud, fo=l. The esti- 
mation of the filling factor commonly uses a comparison 
between the beam solid angle (rts ) with the source solid 
angle (f^c), such that fc — i^c/^B, but other geometric 
factors can be included. 

Since the HI cloud shows a velocity gradient across the 
region (see Figure 2), we will assume that at the velocity 
of the absorption feature (Vlsr — H km s^^), we only 
have to take into account the effect of the absorbing gas, 
neglecting the contribution to the emission (iniTex ^0), 
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since this occurs mostly at other velocities. Then, the 
transfer equation at the center of the absorbing line is, 

TL(0) = (-/cr,)(l-e-^-(°)). 

The resulting optical depth value is obtained using the 
expression: 

1^^(0)1' 



MO) 



In 1 



fcT, 



The beam-averaged brightness temperature of the con- 
tinuum can be computed from the total flux density 
(Sc —72 mJy), and the deconvolved angular size mea- 
sured in the 1.4 GHz continuum image {9c — 3''5), by 
the expression T^ = 688(5^/^7^), with Qc = lAWl, ob- 
taining Tc=3,600 K. The filling factor is f^ ~ ff^/el, 
where 6% is the beam size (-- 15''4), thus fc=0.05. The 
observed line brightness temperature, at the center of the 
absorbing line Tl(0), for Sl(0)=-69 mJy and 6'b=15"4, 
is —177 K. Then the mean optical depth has a very large 
value tl{0) —5. As we note, the derived values for Ti(0) 
and fcTc are similar, such that the ratio Ti(0)/fcTc will 
be very close to unity, making the estimation of the op- 
tical depth very uncertain. Moreover, If we include the 
uncertainty in the estimated angular size and tempera- 
tures, the error in the mean optical depth will be around 
60%, Ti(0) ~5±3. 

A lower limit to the HI mean density can be set by 
assuming an optical depth, tl{0) >> 1. This assump- 
tion is supported by our previous estimate, tl{0) '^ 
5, but with large uncertainties. Then taking simply 
tl(0) >1, Tcx ^300 K and the FWHM of the line 
AV=2 km s~^ (see Table 1) we can estimate a lower 
limit for t he HI column density, computed in the stan- 
dard way (iRohlfsfcWilsonI [20031) . 
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givmg ^Hi >l.lxlO^^ cm^^. This column den- 
sity is consistent with t hat expected for PDRs 
(jHoUenbach fc TielensI Il997f) . We will further assume 
that the physical depth of the region producing the HI 
absorption is approximately one half of the diameter of 
the HI in emission (^42"; ~ 0.2 pc). Therefore an esti- 
mate of the average neutral hydrogen density, < nHi >— 
N/f//L, where L~0.1 pc, is < uhi > >3.5xl0^ cm~^. 
The neutral hydrogen mass associated to the HI emis- 
sion can be obtained assuming a spherical cloud with 
an HI radius of 21 (0.1 pc) and average hydrogen den- 
sity < uhi > >3.5xl0'^ cm"'^, from the expression 
M{HI) — ^R^ < riHi > m.Hi, where R is the radius of 
the HI emission region, and mni is the hydrogen atomic 
mass. We obtain that the HI mass M(HI)>0.4 Mq. This 
value should be taken as a lower limit, si nce the aver- 
age hy drogen density is also a lower limit. iGomez et al.l 
()1998() used a line emission model for the PDR to es- 
timate the total HI mass, obtaining a value as high as 
5 M0. Molecular mass estimates, derived from CO, indi- 
cate that GGD 1 4 has a molecular mass in the range from 
640 to 1500 Mq (iRidge et al.ll200llHiguchi et alJl2009l) . 
which means that the PDR mass is very low (<1% of the 
total gas mass). Of course, the PDR around G213.880- 
11.837 extends only over ~1 , while the CO is detected 
over ~5' (|Higuchi et al.ll2009f ). 



4.2. Kinematics of the HI gas in GGD 14 

The high optical depth toward the CHII region 
(■'"l(O) >1) supports the existence of high density atomic 
gas, where the cometary CHII region (G213. 880-11. 837) 
is embedded. The spatial extension of the cometary 
HI 21-cm emission (~42 ) is ~10 times larger than 
the ionized CHII region (<3.5" at 1.4 GHz). In the 
champagne model a strong density gradient is required, 
with the head of the cometary CHII region embedded 
in the densest part of the molecular clump and the 
ionized gas expanding a s ymm etrically out of the dense 
clump (JGarav fc Lizanol Il999l). For G213. 880-11.837, 
iGomez et al.l (|1998[ ) interpreted the kinematics of the 
ionized gas, derived from the H92a line, as a cham- 
pagne flow due to the presence of a velocity gradient mea- 
sured from the head (located to the NW with Ylsr '^H 
km s~^) to the tail (located to the SE with Vlsr ^^15 
km s^^). 

Recent numerical results, show that an expanding H II 
region should be an efficient trigger for star formation 
in molecular clouds if the mass of the ambient molec- 
ular material is large enough (JHosokawa fc Inutsukal 
|2006[ ). In the particular case of the GGD 14 com- 
plex, a dense molecular region should exist in the 
vicinity of G213. 880-11. 837, where recently a cluster 
of low-mass pre-main-sequence stars has been found. 
In addition to the radio continuum co mpact sources 
(jGomez. Rodriguez fc Garavll2000l I2002D. there are in- 
frared studies in the 2-100 /im ran ge (jHarvev et al"]|1985l : 
lHodappl[l99l iFang fc Yaol 120041 ) that provide evidence 
for a recently fori ned cluster of sourc es in the vicinity of 
G213.880-11.837. iFang fc Yaol (|2004D estimated the age 
of this cluster in the order of 2—5 x 10^ yrs. 

The molecular gas detected toward GGD 14 shows a 
complex velocity structure with the presence of a CO 
molecular outflow with the major axis along the NW- 
SE direction (P. A. =—60°), and centered near the ra- 
dio source VLA 7. In particular, we note that, i n 
the C0(2^.1) images (|Little et a]||l990HOin et al.ll2008D . 
there is strong self-absorption of this line at Ylsr "^H 
km s~^, which c orrespond to the ambient molecular gas 
velocity. iQin e t al. (2008) shows a 1^00(2^1) spec- 
trum toward this region (not affected by self-absorption), 
that peaks at a velocity of ~11.5 km s~^. In what fol- 
lows we will assume this value for the am bient molecular 
gas toward G GD 14. Th ere is also NH3 ('Torrelles et alj 
il983; Glisten fc Marcaid e 1986; Anglada ct al. 1996), CS 
(lAnglada et aLlll996l ) and HCO+ emission (JHeaton et aTl 
Il988f ). toward GGD 14, with velocities in the range from 
9 to 12 km s~^. These results point out that there is a 
dense region in the vicinity of the CHII region that could 
be responsible of the champagne morphology observed in 
both the CHII and the HI region. 



4.2.1. Champagne HI flow around G213. 880-11. 837? 

Figure 3 shows a near-IR continuum image at 4.5 /im of 
Spitzer Space Telescope^ toward GGD 14. We used the 
processed Spitzer IRAC (jFazio et al.l l2004l ) Basic Cali- 
brated Data (PID: 6. PI: Giovanni Fazio), which are 

■^ The Spitzer Space Telescope, is operated by the Jet Propulsion 
Laboratory, California Institute of Technology under a contract 
with NASA. 
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available in the archive of the Spitzer Science Center. In 
this image the stellar cluster is clearly appreciated. We 
overlay in this image the with integrated HI 21 cm line, 
from 1 to 20 km s~^. The stellar cluster is mainly lo- 
cated toward the HI absorption region. We notice that 
the compact radio sources VLA 1 (G213. 880-11. 837) and 
VLA 7 h ave strong 4.5 iim. counterparts no t observed at 

2.2 fj,Ta (|G6mez. Rodriguez fc Garavl |2002|) . which sug- 
gests that these sources are deeply embedded. The 4.5 
fim counterpart toward VLA 1 is first reported in this 
work, and that t oward VLA 7 was recently detected by 
iSatoet al.l(l2008D. T he near-IR source IRS 9Mc reported 
bv lSato et al.l (|2008[ ) is localized near VLA 7, and shows 
a bipolar infrared nebula aligned in the same direction as 
the bipolar molecular outflow (NW-SE; 'Rodrigue z et al.l 
[1982:) . We also noti ce that several of the compact radi o 
sources reported bv lGomez. Rodriguez fc Garavl ()2002D . 
indicated with crosses, have 4.5 /im counterparts. 

Figure 3 shows the velocity integrated HI 21 cm line 
toward G213. 880-11. 837, which exhibits a cometary like 
morphology, resembling (at much la rger scales) that of 
the ionized gas at high frequencies (iTofani et al] 119951 : 
iKurtz et al.lll994HG6mez et al.lll998[) . The limb bright- 
ened PDR expected at the head of the cometary re- 
gion G213. 880-11. 837 is not observed. This can be due 
to the HI absorption feature that complicates the full 
imaging of the PDR or to the low angular resolution 
of the HI observations. In the latter case, one can set 
a lower limit to the HI volume density at the head of 
the cometary. Assuming a PDR size equal to the beam 
size, I = 1.5 X 1O*AU(0b/15 ){d/lkpc), and a column 
density of the PDR Nhi -^ 3 x 10^^ cm'^ such that 
the dust opacity is of order unity, nni > Nhi/1 — 

1.3 X 10^ cm""^, which is a reasonable density for molecu- 
lar clouds. Molecular line observations with high spatial 
resolution are required to try to detect this high density 
molecular gas at the head of the cometary. 

In order to assess whether or not the HI emission is pro- 
duced by a champagne flow, we made position-velocity 
diagrams of the HI emission along two perpendicular di- 
rections, centered on the CHII region (peak continuum 
emission). The cuts are at position angles of -1-25° and 
—65° which correspond to the orientations of the mi- 
nor and major axis of the ionized cometary CHII region 
([Gomez et al.|[l998l) . Figure 4 shows the PV diagrams of 
the HI gas along these two directions. They both exhibit 
absorption and emission features. The spatially unre- 
solved HI line in absorption seen at the position of the 
CHII region is strikingly broad in velocity, ranging from 
5 to 13 km s~^. We also see that at this position the 
HI emission reaches its maximu m redshifted velo cities 
with a value of about 18 km s~^. iQin et all (|2008f) have 
shown that the velocity of the ambient molecular gas to- 
ward GGD 14 is ^^11.5 km s^^. We interpret these two 
features as indicating the presence of an expanding shell 
of HI gas around the CHII region, with an expansion 
velocity of about 6.5 km s^-^. 

The PV diagram along the position of —65° shows that 
the emission is very extended toward the SE direction, 
reaching velocities in the tail of about 14 km s~^. Since 
the ambient velocity of the molecular gas is 11.5 km s~^, 
we interpret the PV structure seen SE of the CHII region 
as indicating a champagne flow of the HI gas. In a cham- 



pagne flow the gas accelerates to 2 — 3 times the sound 
speed (see Table 2 of lShu et al.ll2002D . Since the observed 
velocity shift in the i onized gas is Az; ^ 15 km s~^— 11 
km s~^ = 4 km s~^ ([Gomez et al.lll998D while a shift of 
Av '-^ 20 — 30 km s~^ is expected, this ionized flow can 
be understood as a champagne flow only if the source 
is inclined close to the plane of the sky, with an angle 
i < 1 1°, and has a small opening angle. iGomez et al.l 
(|1998() correctly interpreted the kinematics of the ion- 
ized gas as a champagne flow but they did not estimate 
the inclination angle of the flow with respect to the plane 
of the sky, they only suggest that the tail should be going 
away to explain the redshifted velocities. In the present 
interpretation we propose that in fact the orientation of 
the ionized champagne flow should be close to the plane 
of the sky with the tail having the redshifted velocities. 
On the other hand, in a PDR, one expects gas tempera- 
tures of T ~ 500 K, i.e., sound speeds ~2 km s~-^. Thus, 
for a champagne flow in the HI neutral gas, the expected 
velocity shift is of the order of ~4 — 6 km s~^. Thus, the 
observed velocity shift, Avhi ^ 3 km s~^, is consistent 
with a champagne flow of the HI gas, with a large opening 
angle such that the observed velocity in the line-of-sight 
is small. Note that the scale of this HI flow is ^10 times 
larger than the ionized champagne flow, therefore, the 
collimation of this HI neutral gas is probably much lower 
supporting the presence of a larger opening angle. 

In contrast, in the cut along the perpendicular direc- 
tion to the symmetry axis passing through the HII region 
(P.A.=-f25°), the HI gas in emission is seen at both sides 
of the cometary head and no velocity gradient is appreci- 
ated. High angular resolution molecular observations are 
needed to confirm this kinematics. Another HII region 
that shows a cometary morphology in both the contin- 
uum and HI emission is Gill. 61-1-0. 37, where a cham- 
pagne interpretation was also proposed (Leb ron et al.l 
[200ih . 

Hydrodynamica l and radiative transfer models by 
iHosokawal ([20071 ) show that an inhomogeneous medium 
can change the time evolution of expanding HII regions 
and PDRs (as compared with the evolution in a homo- 
geneous medium). In a medium with a radial density 
profile of the form n{r) oc r^™ with < w < 1, the 
dense shell formed by the shock front shields the far- 
ultraviolet (FUV) radiation and prevents the expansion 
of the dissociation front beyond the shell. The CHII re- 
gion is surrounded by a cold molecular envelope in this 
case. In a stepper density gradient with w > 1, as ex- 
pected in a "champagne flow," shielding of the FUV radi- 
ation becomes inefficient and an extended PDR develops 
around the HII region. However, only ii w > 1.5, the H2 
dissociation front propagates through the ambient cloud. 
Since in G213.880-11.837 and Glll.61+0.37 we observe 
extended PDRs around the HII regions, it is likely that 
steep density gradients exist in both objects. 

5. CONCLUSIONS 

We present high angular resolution (15 ) VLA obser- 
vations of the HI 21 cm line toward the compact HII 
region G213.880-11.837. The HI gas, produced by a pho- 
todissociated region around the CHII region, exhibits a 
cometary morphology with an extension, along the major 
axis, of ^45 . The kinematics and morphology of the HI 
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gas around G213. 880-11. 837 indicates that the neutral 
hydrogen gas is in expansion in a similar way as the ion- 
ized gas. In both cases the gas accelerates toward lower 
density regions and redshifted velocities (from ~11 to 14 
km s~^). The peak position of the HI absorption coin- 
cides with the peak continuum emission of the CHII re- 
gion (VLA 1; G213. 880-11.837). This is the second (after 
Gill. 61-1-0. 37) cometary HII/HI region known. These 
observations suggest the existence of champagne flows in 



PDRs associated with cometary HII regions. 
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Fig. 1. — Integrated spectrum of the HI 21 cm line emission toward GGD 14. 
Gaussians (see Table 1). 
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Fig. 2. — Individual channel images of the 21 cm HI emission and absorption toward the GGD 14 region. One out of every two channels 
was plotted. Contour levels are -20, -15, -10, -7, -3, 3, 5, 7, 10, and 15 X 2 mjy beam~^, the rms noise of the images. The continuous and 
dashed contours indicate emission and absorption, re spectively. The cross and the trian gle mark the position of the radio continuum peak 
at 1.4 GHz of VLA 1 (G213. 880-11.837) and VLA 7 l IGomez. Rodriguez fc Garavll2000l 'l. respectively. 
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Fig. 3. — Superposition of the near-IR Spitzcr image at 4.5 fj,ra (gray scale) on the VLA integrated HI hne (from 2 to 20 km s ^) with 
solid (emission) and dashed (absorption) contours (-20, -10, -3, 3, 6 and 9 X 0.7 mjy beam~^ km s"'^). The large X symbol indicates the 
radio continuum peak position of the HII region (VLA 1) and the triangle the peak position of VLA 7. The crosses mark the peak position 
of the compact radio sources reported by IGomez. Rodriguez fc Garay 1 2002ll . 



An expanding HI PDR flow in GGD 14 




LSR VELOCITY (km/s) 



Fig. 4. — Position velocity diagrams of the HI line toward GGD 14. Top: cut along the symmetry axis of the cometary HI structure 
(P.A.=— 65°). Bottom: cut along a perpendicular direction to the symmetry axis passing through the HII region (P.A.=+25°). Contour 
levels are -24, -21, -18, -15, -12,-9, -6, -3, 3, 6, 9, 12, 15, and 18 X 2.5 mjy beam~^, the rms noise in these images. The insets shovir the 
direction of the cuts writh respect to the cometary morphology. 
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TABLE 1 

GAUSSIAN COMPONENTS 

IN THE FIT TO THE HI 

LINE PROFILE 



Si Vlsr AV_ 

(injy) (km s"-*-) (km s^-*-) 

-98±36 1.2±0.6 3±1 

-49±24 6.4±1.5 5±3 

-69±27 11.1±0.6 2±1 

267±18 15.5±0.1 4±1 



